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ABSTRACT 

We have obtained a series of high-resolution optical spectra for 2MASSW 
J1207334-393254 (2M1207). Two consecutive observing nights at the ESO Very 
Large Telescope with the UVES spectrograph yielded a timeseries with a res- 
olution of ~ 12min. Additional high-resolution optical spectra were obtained 
months apart at the Magellan Clay telescope using the MIKE instrument. Com- 
bined with previously published results, these data allow us to investigate changes 
in the emission line spectrum of 2M 1207 on timescales of hours to years. Most of 
the emission line profiles of 2M 1207 are broad, in particular that of Ha, indicating 
that the dominant fraction of the emission must be attributed to disk accretion 
rather than to magnetic activity. From the Ha 10 % width we deduce a relatively 
stable accretion rate between 10^^^'^'"^^'^Mq/jt for two nights of consecutive 
observations. Therefore, either the accretion stream is nearly homogeneous over 
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(sub-) stellar longitude or the system is seen face-on. Small but significant vari- 
ations are evident throughout our near-continuous observation, and they reach 
a maximum after ~ 8 h, roughly the timescale on which maximum variability is 
expected across the rotation cycle. Together with past measurements, we confirm 
that the accretion rate of 2M 1207 varies by more than one order of magnitude 
on timescales of months to years. Such variable mass accretion yields a plausible 
explanation for the observed spread in the M — M diagram. The magnetic field 
required to drive the funnel fiow is on the order of a few hundred G. Despite 
the obvious presence of a magnetic field, no radio nor X-ray emission has been 
reported for 2M 1207. Possibly strong accretion suppresses magnetic activity in 
brown dwarfs, similar to the findings for higher mass T Tauri stars. 

Subject headings: Accretion, Accretion Disks - Stars: Low-Mass, Brown Dwarfs 
- Stars: Circumstellar Matter - Line: Formation, Profiles - Stars: Planetary 
Systems 



Introduction 



In recent years, there is increasing evidence that young brown dwarfs evolve similarly to 
their higher- mass counterparts, the T Tauri stars (TTS). In particular, their formation pro- 
cess includes a phase of mass accretion from a circums tellar disk, whose presence is probed by 



infrared emission in excess of the photosphere (e.g., iNatta et al.ll2002l : IJayawardhana et al. 



2003al : iLuhman et al.ll2005l : lAUers et al.l 120061 ). Spectroscopic features, such as the strength 
and shape of Ha emission, directly trace the accreting material. In the stellar regime, 'weak 
TTS' have no optically thick circumstellar matter, and their Ha emission is purely chromo- 
spheric. In 'classical TTS', the line emission is dominated by contributions from accretion. 
Next to displaying much stronger lines, thei r line profiles are broad er due to the high velocities 
in the accreting gas and Stark broadening (IMuzeroUe et al.ll200ll ). Winds may play a role as 
well, and are recognized by bluesh i fted absorption components and for bidden emission lines 
(jMundtlll984l : Edwards et al.lll987l : lHamanrull994j : iHartigan et al.lll995l ) . Analogous observa- 
tions have been made for brow n dwarfs, where classi cal TTS-like Ha emission persists down 
to the deuterium burning limit (IMohanty et al.ll2005l ). Forbidden emi ssion lines, tracing out- 
fiows have been reported from a small number of brown dwarfs (e.g.. [Fernandez fc Comeron 
200ll : iBarrado v Navascues et"aDl2004j l 



Brown dwarfs are intrinsically faint objects, therefore a detailed study of magnetic activ- 
ity and accretion in young brown dwarfs is restricted to the most nearby ones. A particularly 
well-suited playground is the TW Hya association (TWA), with four confirmed brown dwarfs 
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at an age of ~ 8 Myr and a distance of ~ 50 pc. 2M 1207 is the best-studied among the brown 
dwarfs in the TWA. Its youth was confirmed by low gravity and Lithium absorption, and its 



TWA membership by rad ia l velocity a r id pro per motion studies (lGizisll2002l : iMohanty et al. 



20031 : IScholz et al.ll2005bl ). iMamajekl (120051 ) has revised the previous distance estimate of 
2M 1207 using the moving cluster method, resulting in 53 ± 6 pc. 2M 1207 is the oldest 
brown dwarf known to actively accrete from a disk, thus representing a benchmark in con- 
straining the disk lifetimes in the substellar regime. The object has recently aroused stron g 
interest because of the discovery of a planetary mass companion (jChauvin et al.ll2004l . |2005| ). 
This companion is at a separation of 40 AU and not expected to influence the magnetic and 
accretion activity of the brown dwarf. 

Stellar astronomers are intrigued by 2M 1207 being the closest known brow n dwarf with 
a disk. No substantial L' excess was detected by IJayawardhana et al.l (l2003ar). bu t mid-IR 
photometry established the presence of circumstellar material (jSterzik et al.ll2004j ). This is 
also corroborated by Spitzer photom etry, where 2M 1 207 is detected in all IRAC bands and 
at 24 /im above photospheric levels (IRiaz et al.ll2006l ). The spectral energy distribution of 
2M 1207 shows evidence for dust evolution, similar to the higher-mass disk-bearing TTS in 
the TWA. 2M 1207 is therefore a lower-mass analog of the classical TTS. 

Strong Ha e mission ha s been reported right with the identification of 2M 1207 as a 
TWA candidate (lGizisll2002l ). Subsequent high-resolution spectra have shown that the Ha 
profile is broadened a nd asymmetric, thus displaying characteri stic signatures o f accr etion 
( IMohanty et al.ll2003l ). The Ha equivalent width me asured by iMohanty et al.l (120031 ) was 
only about one-tenth of the earlier value obtained by Gizis ( 20021 ). and this variability was 
interpreted as a sign for variable mass accretion rate. In a subsequent study of emission line 
variability, aiming at constraining these variations, 2M 1207 presented dramatic changes i n 
the shape and intensity of its Ha line in data obtained two months apart (jScholz et al. I l2005ah . 
The observed profile variability was consisten t with a scenario of asyrnmetri c accretion with a 
direct view into the (magnetic) funnel flow (jScholz fc Jayawardhana 2006 ). In this picture, 
variability is also expected on the time scale of the rotation period (~ Id), as the hot 
accretion spots are corotating with the star. Existing data have remained inconclusive in 
this respect due to poor sampling. Therefore, we have collected new high spectral resolution 
data with the aim to monitor the Ha variability on shorter time scales. 

Here we present the results obtained during two consecutive nights of observation, and 
from five additional spectra taken weeks to months apart. The observations are described 
in Sect. HI In Sect. Owe present a detailed analysis of the characteristics of the Ha emission 
and other prominent emission lines in the optical spectrum of 2M 1207. The results are 
discussed in Sect. HJ where the spectral variability is examined on various timescales from 
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hours to years. Estimates for the accretion rate and the magnetic field strength of 2M 1207 
are presented, and put in context with other (sub) stellar objects. Finally, we compare the 
accretion and activity characteristics of all known brown dwarfs in TWA. Concluding remarks 
are found in Sect. O 



2. Observations 

2M 1207 was observed with the Ultraviolet and Visual Echelle Spectrograph (UVES) at 
the ESO VLT2 during two nights May 8-9, 2006 in the dichroic mode that covers the full 
spectral range of UVES by simultaneous use of both the blue and the red arm. We used the 
standard setting DICHR7^2 centered at 437 nm for the blue arm and at 760 nm for the red 
arm. This setup yields simultaneous spectral coverage from 3730 — 4990 A on the blue chip 
and from 5650 — 9460 A on the red mosaic. A slit width of 0.8" was chosen resulting in a 
resolving power R ^ 50000. The CCDs were rebinned to 2 x 2 pixels. 

The integration time was 600 s. During the two observing nights, we obtained a total 
of 59 exposures of 2M 1207 at blue and red wavelengths. Except for some constraints due to 
calibrations and weather the frames were taken consecutively, such that the time resolution 
for the investigation of spectral changes is ~ 12min over a large part of the observations. 

T he data were process ed with the UVES pipeline implemented in the MIDAS environ- 



ment (IBallester et al.l 120001 ). Reduction steps included bias subtraction, fiatfield correction, 
and wavelength calibration. The final one-dimensional spectra have been corrected for con- 
taminating sky light and cosmic rays have b een removed. Each spectrum was corrected 



for the radial velocity of 2M 1207 (11.2 km/s; iMohanty et al.ll2003r ) and for the barycenter 
motion of the Earth. 

We extracted a spectrum from each of the individual science frames. In addition, to 
increase the S/N for the detection and measurements of weak spectral features, we built 
averages of each four consecutive frames. Thus, the effective exposure times of these merged 
spectra are 40min. Three individual frames with poor S/N were not considered in the 
merging process, such that we obtain 14 average spectra from the two nights. The spectra 
were rebinned to R = 50000 around individual emission lines before the scientific analysis 
was carried out. 

In February and April 2006, we obtained five additional exposures of 2M 1207 using the 
Magellan Inamori Kyocera Echelle (MIKE) spectrograph at the Magellan Clay telescope on 
Las Campanas, Chile. These data were reduced using standard routines within IRAF, in- 
cluding bias subtraction, fiatfielding, wavelength and relative fiux calibration, and removal of 



- 5 - 



cosmic rays. For more details of the reduction procedure, see lScholz &: Jayawardhanal (120061 ) . 
Exposure times for the MIKE spectra were typically lOmin, and the nominal resolution is 
R ^ 25000. 



3. Emission line spectrum 
3.1. Ha emission 

The main aim of this study is to examine the variability of the Ha emission of 2M 1207 
with good time resolution and on time-scales of ^ 1 d, corresponding to its suspected rotation 



period. An estimate for the period is given bvlScholz et al 
velocity of 2M 1207 is v sini = 13 ± 2 km/s (IMohantv et al. 



(l2005al) . In brief, the rotational 



20031) . Together with a radius of 



0.27 -Rq extracted from evolutionary models (IChabrier et al.ll2000l ). the period is expected to 
be P/ sin i ^ 25.2 h. Assuming a random position for the inclination, the most probable value 
for sini is 0.81, and for this va lue the rotation period is 20.4 h. Magnetospheric accretion 
models (IHartmann et al.l Il994l ) have shown that the presence of a red absorption feature 
in the Balmer lines (frequently seen in 2M 1207; see below) indicates high inclination. For 
i > 60° the expected period is > 22 h. 

The Ha emission of 2M 1207 is thought to be dominated by accretion: the broad wings 
with velocities of up to ±200 km/s can not be explained by magnetic activity, the rotational 
velocity is comparatively small and can not be responsible for the br oad wings, and signatures 
for winds are very weak (IMohanty et al.l l2005l : IWhelan et al.l 120071 ) . 



The most widely used indicator for identifying accreting stars is the equivalent width 
Weq of the Ha line, where weak emission is due to chromospheric magnetic activity and 
strong emission, above a certain threshold of Weq, is attributed to accretion. However, 
its dependence on the continuum level precludes Weq as a universal measure. To take 
account of the declining continuum for later spectral types, the threshold for Weq that 
separates accretors from non-accretors among TTS is usually adapted in discrete steps from 
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Furthermore, equivalent width measurements of faint continuum sources 
such as brown dwarfs are intrinsically associated with considerable uncertainties. This is 
due to the low S/N and the presence of overlapping absorption features that form a 'pseudo- 
continuum' w hich depends on the s pectral resolution. To avoid these problems, following a 



suggestion bv lWhite fc Basril (120031 ) . the full width at 10 % of the peak height {Wio%) of Ha 



has been established as more reliable accretion indi cator: Objects with Wiq% > 200 A can 



be considered accretors (jjayawardhana et al.ll2003bl ) 
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3.1.1. Profiles and line widths 

We use lVio% as the primary accretion indicator, but examine also other parameters that 
characterize the hne intensity and shape. Fig. [1] shows the time-series of these parameters 
obtained during our run with UVES on May 8-9, 2006. The 10 % width is clearly above the 
accretion threshold throughout our observation (H^io% ~ 280...320 km/s); on average it is 
higher during the second night, in the course of which it slightly declines. The equivalent 
width has ups and downs during the first night, and has its maximum at the beginning of 
the second night, after which it declines dramatically and in a continuous way until the end 
of the observation. The errors in VFeq are dominated by the (faint) continuum level. We 
estimated the uncertainty of the continuum flux by varying the sky extraction area, and 
found that it is of the order of 15%. Then we computed error bars for PVeq assuming for 
the continuum flux 0.85 ■ /cont and 1.15 ■ /cont, respectively. Fig. [1] shows that when taking 
into account these uncertainties, the observed variations of IVeq remain significant. 

The Ha emission is double-peaked throughout all our observations from February to 
May 2006. Fig. [2] displays the time-series of all 14 normalized average line profiles /n(A) 
from the UVES run (May 2006), and Fig. [3] shows the 5 normalized profiles obtained with 
MIKE (February and April 2006). For clarity, in Figs. [2] and [3] the individual spectra have 
been shifted with respect to each other in the vertical direction. 



Scholz et al.l (l2005a] ) have shown that the double-peaked Ha profile of 2M 1207 can 
be interpreted as a broad emission line onto which a red-shifted absorption component is 
superposed. For the minimum of the absorption feature we measure a position that varies 
in the range Afmin ~ —10... + 70 km/s with respect to the expected wavelength (Afmin ^ 
+30... + 70 km/s in February 2006, ^v^^^ +10 km/s in April 2006, and l^v^^^ ^ -10... + 
40 km/s in May 2006). From Fig. [1] it appears that twice during the extensive monitoring 
in May 2006 fmin undergoes a discrete jump from high redshift towards the expected line 
center. These times correspond to line profiles in which the absorption trough is wide and 
asymmetric, such that the position of the minimum can not be determined precisely (Fig. [2]). 
As a measure for the asymmetry of the emission component. Fig. [T] shows the difference 
between the 10 % width on the blue (left) and on the red (right) side of Ai; = 0. In the 
course of the UVES observations, Vrio%(^ — ^) varied from —40... + 40 km/s. In particular, 
during the second night the line moved systematically from the blue to the red. This is in 
contrast to the absorption feature shown above to be redwards of Aw = throughout the 
observation (cf. also Fig. [2]). 

We examined the shape of the Ha line and its change in time by modeling each of the 
14 average profiles with two quasi-Gaussians, one with positive normalization representing 
the emission component and one with negative normalization for the absorption component. 
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Each of these has three free parameters (position, width, and normahzation). The equivalent 
widths of the emission WEQ-em and the absorption WEQ;abs component derived from the fits 
turn out to be strongly correlated (Fig. H]). Since WEQ;em is by definition 'normalized' to 
the continuum this relation is not trivial. In Fig. H] we also plot the observed flux, i.e. the 
difference between the flux in the emission profile and the absorption profile, versus the 10 % 
width. 

The Ha 10 % width of 2M 1207 has been measured at eight epochs so far (see Table [T]). 
In most of the high- resolution spectra collected so far 2M 1207 showed Wio% clearly above the 
accretion threshold, although in some occasions there were only marginal signs of accretion 
(1^10% ~ 200km/s). 



3.1.2. Quantifying variability 

To quantify the Ha variability as a function of wavelength we examined the normalized 
variance profile 

^n(A) = a^(A)//„(A), (1) 

where the varia nce profile a'^{\) = Xlili [-^n,i(A) — /n(A)]^ measures flux variations across 
the line profile ( Johns &: Basri 19951 ). In the bottom panels of Figj2]we plot a/o^ for three 



distinct parts of the observation in May 2006. From left to right, are shown the variances 
obtained from all 19 spectra from the first night, the first 20 and the last 20 spectra of 
the second night, respectively. This way, the individual v^of give information about the 
variations as a function of wavelength on short time-scales (within the ~ 4 h covered by 19 
respectively 20 consecutive frames), and a comparison of the different -y/of's probes changes 
up to ~ 1 d. 

The normalized variance in the continuum expected from Poisson noise is cr^exp(A) = 

o"„(A)(\/7//)^. The zero-variability level \/o^^^ is shown as dashed line in the bottom 
panels of Fig. [21 Its comparison with the \/af 's shows clear signs for variability in the line 
profile. Throughout the time series, the variance profile (i.e. variability as a function of 
wavelength) features three distinct peaks, roughly at positions of — 100... — 150, —20... + 20, 



and -t- 100... + 150km/s. We note that the same structure was seen by lScholz fc Jayawardhana 



(120061 ) in their Ha time series of 2M 1207. The position of these peaks (indicating the 
velocities of maximum changes in the Ha flux) is roughly coinciding with the position of the 
two maxima and the minimum in the Ha profile. 



The variance profile yields only the time average of the changes as a function of wave- 
length. To study the profile variations in more detail we examine the change of the flux 
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between consecutive spectra across the line. The relative flux change is given by 

where /n(-^) is the normalized profile, z is a running number denoting the exposures and 
j = i + 1. These relative flux changes p(X) are used to inv estigate variations as a function 



of the time lag between two spectra (j Johns fc Basril Il995l ). To this end, the i5(A)'s are 



computed not only for two subsequent spectra but for all pairs of the 59 spectra from May 
8-9, 2006, i. e. for each time lag 6t = t2 — tl > sampled by our UVES data. Then, we 
group the D(A)'s in time lag bins of Ih (At) and calculate the standard deviation crAt(A) of 
the D(A)'s in each of these bins. Finally, we calculate for each time bin the average of these 
standard deviations over the whole line profile (o"At)- The result is shown in Fig. El where 
the numbers on top of the panel indicate the number of data points in each of the time lag 
bins. The gap is produced by the absence of observations with time lags between ^ 9...16h, 
i.e. the data to the left of this gap corresponds to spectra from the same night, and the data 
to the right corresponds to spectra from two different nights. The level of profile variability 
doubles from ~ 2% for the shortest time lag (Ih) to a maximum of ~ 4%, that is reached 
for a timelag of ~ 8 h. 

To test whether the variations are correlated across the line profile we computed the 
correlation coefficient for each pair of spectral bins using the profiles ^^(A), with i = 1, ...,N. 
This results in a correlation matrix ri j of dimension A^b x where A^b is the number of 
spectral bins. The confidence level of each coefficient in the matrix is computed using the 
error function 

p _ I n,j I ViVb . . 

such that 1 — -Eij is the confidence level for the correlation to be true. 

The correlation matrix has been computed for different subgroups of Ha profiles from 
the UVES run. Fig. [6] displays the most interesting correlation matrices j as contour plots. 
Only values with confidence level > 99.9% are considered in this representation. The first 
two diagrams from the left represent the matrices of each 19 consecutive spectra, obtained 
in the first and the second observing night, respectively. The cross-shaped pattern seen in 
the leftmost diagram is peculiar. The center of the cross is offset by ~ +25km/s to the 
red side, and it represents the spectral region near the minimum of the line profile (see 
Fig. [1]). Closer examination shows that this pattern stems from a short time interval. This 
is demonstrated in the rightmost diagram of Fig. [6], that displays the correlation matrix for 
exposures i = 7, 10, roughly corresponding to the second average profile in the leftmost 
panel of Fig. [2l Clearly, during this time interval, the profile changes near the absorption 
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minimum are uncorrelated with all velocities except for the velocities in close vicinity. The 
double-square sh aped pattern observ ed during the second night (middle panel of Fig. [6]) is 
more typical (cf. iJohns fc Basrilll995l ). The absence of contours in the upper left and in the 
lower right of this graph means that red and blue wings vary in an uncorrelated fashion. 
The range near v ~ -|-75km/s, corresponding to the position of the right emission peak, is 
particularly well correlated with the whole line. 



3.2. Other signatures of accretion 

Our high S/N UVES spectra allow us to present the first time series for other optical 
emission lines of 2M 1207. In particular, we detect the Balmer lines up to He. This latter 
one is partly blended with CaIIHA3968. 

The shape of the H/? profile follows closely that of Ha discussed in detail above. All 
higher Balmer lines are predominantly characterized by a double-peaked profile, but at 
times turn into a single-peaked shape. This is demonstrated in Fig. [TJ that shows the 
contemporaneous Balmer series for two epochs. Each panel in Fig. [7| is obtained from an 8- 
frame average. When present, the absorption minimum is found at roughly the same position 
for all Balmer lines, but the ratio of the emission fluxes left and right of the minimum 
decreases systematically for the higher Balmer lines, i.e. the lower Balmer lines are red- 
dominated and the higher Balmer lines are blue-dominated. 

In contrast to the Balmer lines, the Ca II H+K emission is narrow with FWHM of only 
~ 20km/s (see Table [2]). Furthermore, we detect several lines of He I (A4471, A5876, A6678, 
A7066) and the Call infrared triplet (hereafter CaIRT) in emission. All these lines are 
predominantly detected in accretors and, therefore, thought to be related to the accretion 
process. The Ca II A8 662 flux was shown to be directly related to the mass accretion rate 



( Mohanty et al.ll2005l ). The three lines of the Ca IRT appear with roughly similar strength 



inconsistent with formation in an outflow (IReipurth et al.lll986l : [Fernandez &: Comer6rul200ll ). 



No Na D emission is see n, consistent with the prediction of magnetospheric accretion models 



for low accretion rates (IMuzeroUe et al.ll200ll ). 



We measured H^eq and the FWHM of the stronger ones of the emission lines discussed 
above. The continuum was estimated from two 5 A wide line-free regions left and right of 
the respective line. The observed range of VFeq is reported for each line in Table [2] together 
with its approximate uncertainty. Given the often loose use of equivalent widths in the 
literature, it seems appropriate to stress that Weq is associated with considerable errors in 
objects with faint continuum. This is clear from the numbers given in Table [21 The errors 
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cited in that table do not include systematic uncertainties, again affecting mainly the (faint) 
continuum, that we discovered in the process of defining the sky extraction area. Since we 
are interested in relative changes within the UVES dataset, systematic uncertainties are not 
a big concern. However, they should be kept in mind when comparing data obtained with 
different instruments, and analysed in different ways by different astronomers. 

All emission lines are stronger in the second night of the UVES observations, and they 
are particularly weak towards the end of the first night (around MJD 53863.25 = UT 05 h 
on May 8, 2006). At the same moment the shape of some of the strongest emission lines 
underwent a change: (i) He I A5876 temporarily becomes flat-topped and symmetrical, while 
it is markedly asymmetric with an extended red wing throughout most of the two nights 
(see Fig. [H]); (ii) in Ha the 10% width is at its minimum, and the deficiency in the red vs. 
the blue peak is particularly strong; (iii) in the higher Balmer lines the red portion of the 
line disappears almost completely, and (iv) some of the weaker lines nearly disappear at the 
same time. 



4. Discussion 

4.1. Spectral variability 

We have examined the optical emission line variability of 2M 1207 at high spectral 
resolution in data obtained at various epochs through February, April, and May 2006. In 
particular, we present for the first time a near continuous spectral sequence from two consec- 
utive nights with a time-resolution of ~ 12 minutes. The study presented here complements 
previous work that was characterized by much poorer temporal sampling. 



4.I.I. Short timescales (days to hours) 



Scholz et al.l (l2005al ) have attributed the drastic changes of the Ha line profile seen in 



earlier high-resolution spectra of 2M 1207 to varying viewing angle of the accretion column(s) 
over the course of the brown dwarf's rotation cycle. The absorption component seemed to 
appear and reappear on a time-scale of ~ 1 day, an interval approximately coincident with the 
expected rotation period, therefore suggesting that the accretion spots moved with respect 
to the observer. However, with a maximum of 3 spectra per night, the available data did 
not provide adequate sampling across the presumed rotation cycle. 



Within the two nights of near-continuous observation in May 2006 there is clear vari- 
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ability present in the Ha emission and other emission hnes, albeit much less pronounced 
and systematic than expected from previous data. Here, we summarize its characteristics: 
(i) The line profile is double-peaked throughout the two nights, displaying the central ab- 
sorption reversal characteristic for accretion in star-disk systems seen near edge-on; (ii) the 
absorption feature moves between —10... -|- 40km/s; (iii) the absorption minimum widens 
twice during the observations with a time interval of ~ 23 h; (iv) variability across the line 
is inferred from the variance profile and the correlation matrix; and (v) the line flux varies 
by a factor of ~ 6, correlated with the (small) variations of the 10 % width. 

The variable redshift and width of the absorption minimum implies that the viewing 
angle of the absorbi ng layers of the accretion column changed throughout the observation. 



Bouvier et al.l (120031 ) explained the radial velocities of the redshifted absorption components 
of the TTS AATau with changes in the structure of the magnetosphere. In this scenario, 
during times of high observed absorption redshift an observer at high inclination looks flat 
onto the magnetic fleld lines that carry the absorbing material, while low absorption redshift 
implies a large angle between the absorbing material and the line-of-sight. The different 
viewing angles can be realized by inflating or deflating the star-disk fleld. In the case of 
2M 1207, the observed changes that occur on time-scales shorter than the rotation period 
may suggest that the magnetosphere is somewhat inhomogeneous across stellar longitude. 
However, the changes are not as dramatic as seen before. The range of observed redshifts 
in the absorption minimum of 2M 1207 is comparable to that of AA Tau but there is no 
blueshifted absorption. Therefore, for 2M 1207 outflows do not play a signiflcant role in the 
formation of the Ha line. 

Although there is no obvious periodicity, the rotation of the star is the most likely cause 
of the short-term changes of the Ha proflle summarized above. We have shown that the 
fractional proflle changes reach a maximum after a time-lag of ~ 8 h, and they remain high 
on a plateau for larger timelags. (The exact position of the maximum is unclear due to 
the absence of data with timelags between 9 — 15 h.) This timescale is roughly consistent 
with one half the expected rotation period, and therefore with the time interval in which 
maximum changes due to rotational effects are expected. 

If possible contributions from chromospheric activity and winds are neglected, the emis- 
sion line flux is expected to be a measure for the accretion rate, and this is supported by 
its correlation with the 10 % width. This correlation is not trivial because the flt function's 
normalization is independent of its width, unlike in a 'canonical' Gaussian. We conclude, 
from the correlation observed between the equivalent widths of the emission and the absorp- 
tion components (Fig. Hj) that the amount of absorbing material is proportional to the total 
amount of material crossing the line-of-sight. 
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With the exception of the CallH + K hnes, the strong emission hnes (higher Balmer 
hnes, He I and CaIRT) follow the Ha emission closely in their shape and variability, indi- 
cating a common origin, i.e. accretion, and chromospheric contributions are not significant. 
In particular, throughout our observations the lines of the Balmer series are mostly double- 
peaked. Only the higher lines B.S and He at times turn into a single-peaked shape. 



4.1.2. Intermediate timescales (weeks to months) 



The Ha line is characterized by a double-peaked profile in all spectra obtained from 
February to May 2006. In the spectra from February and April the redshift of the absorption 
feature is systematically larger with respect to the observations from May of the same year, 
and consequently the line profiles observed in February and April are also more asymmetric 
than the ones from May. However, in the data f rom 2006 the variatio ns are smoother and 
less dramatic than the ones from 2005, for which IScholz et al.l (j2005al ) reported alternating 
single- and double-peaked profiles within individual observing nights. 

If the previous assumption of a high-inclination system rotating at a rate of ~ 1 d 
is correct, the data from 2006 are compatible with a rather homogeneous distribution of 
accretion elements, i.e. magnetosphere, across the (sub)stellar longitudes, with smaller in- 
homogeneities as discussed in Sect. 14.1.11 In contrast, during 2005 the absorption reversal 
appeared and disappeared within one rotation cycle, implying strong variations in the accre- 
tion structure across the stellar surface. We conclude that the distribution of the accretion 
elements or the structure of the magnetosphere has changed within ~ lyr. In both years, 
2005 and 20 06, the mass accre tion rate, as inferred from the 10% width of the Ha emis- 
sion (see e.g. iNatta et al.l 12004 and Sect. 14. 2p . has shown fluctuations by a factor ~ 2 on 



time-scales of weeks to months (cf. Tabled]). 



4.1.3. Long timescales (years) 



Three years have passed since the first high-resolution spectrum of 2M 1207 was ob- 
tained, and in this period there seems to be a tendency towards en hanced mass accretion 
rate. Note, however, that the large Ha equivalent width measured by iGizid (|2002| ) in a low- 
resolution spectrum counteracts this trend, and fluctuations on shorter timescales (discussed 
above) may be superposed on the long-term evolution. Further observations are required to 
confirm our speculation about the secular evolution of the Ha emission. 



- 13 - 



4.2. Accretion rate and magnetic field 



Natta et al.l (120041 ) have parametrized the accretion rate as a function of W^io% for TTS 
and brown dwarfs. From their Eq. 1, we find an accretion rate of 10~^^'^'"~^'^ Mq/jt for 
the range of Ha 10 % widths measured for 2M 1207 during the UVES run in May 2006. 
In previous measurements, M^^c was as low as ~ lO^^^'"^ Mq /jt, i.e. 2M1207 changes its 
accretion rate by a factor of 2 on the timescale of days, and at least by one order of magnitude 
over the years. 

In Fig. [9] we put 2M 1207 in the context with other accreting (sub)stellar objects. A 
correlation between accretion rate and mass (Mace ~ M°' with a ~ 2) has been reported 
in the literature . In F ig. [9] we display this relation using data from the references given in 
MuzeroUe et al.l (120051 ). The range of Mace inferred from all available data of 2M1207 is 
indicated by large circles connected by a vertical line. 

Various explanations have been proposed f or the correlation between Maec and M, and 



its large (> 2dex) scatter. iPadoan et al.l (120051 ) see in the Mace — M relation a consequence 



of large-scale Bondi-Hoyle accretion. In their scenario, the scatter is produced by the dif- 



ferent physical conditions in different star forming environments. [Alexander fc Armitage 



(120061 ) assumed that the Mace ~ M relation reflects the initial parameters of the disk-star 
system, and they showed that the scatter can be reproduced by a combination of varying 
disk initial parameters and the viscou s evolution of the dis k in the course of which the mass 



accretion rate decreases. According to [Gregory et al.l (120061 ) the Mace — M relation arises nat- 



urally from models of magnetospheric accretion assuming a realistic field geometry. Finally, 



Mohanty et al.l (12005! ) suggested declining disk ionization with decreasing stellar mass as the 



origin of the Maec ~ M relation, implying that m ost of the sca t ter is due to intrinsic variabil- 



ity in accretion rates. This view was disputed by lNatta et al.l (120061 ) who did not find strong 



variations of the near-IR emission lines indicative of accretion in a sample of 14 objects from 



pies of individual obiects, that indeed accretion rates do vary substantially ( 


GuUbrine et al. 


1996; 


Alencar et al. 


2005; 


Scholz et al. 


2005a; 


Scholz & Javawardhana 


2006). We point out 



that the range of mass accretion rates derived for 2M 1207 is comparable to the spread in 
the Maec — M relation. Therefore, variability might indeed explain this spread. Repeated 
measurements of Mace foi' known accretors will be helpful to verify our suggestion. 



An alternative method to derive accretion rates has been proposed by iMohanty et al. 



(120051 ). They established that both classical TTS and brown dwarfs follow a linear correlation 
between the Ca II A8662 emission flux and Mace determined independently from Ha modeling 
or veiling measurements. Following their recipe, we estim ate the continuum flux -Fcont of 
2M1207 at 8662 A from the synthetic DUSTY spectra of lAUard etaP J2000h . We adopt 
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Teff = 2800 K, log (7 = 4.0, and negligible veiling (see discussion in lMohanty et al.ll2005l . for 
justification of these choices). The line fiux is then given by -Fcaii = -f^cont ■ Weq- For the 
highest equivalent width measured in May 2006 this results in l og Fcp u [erg/cm^ / s] ~ 5.1, 
and the fit relation for the low-mass sample of iMohanty et al.l (120051 ) yields an accretion 
rate of Mace ~ 10"^°'^ M0/yr. This is in rough agreement with the value we derived above 
from the Ha 10 % width. We caution, that using the Ca II A8662 line for estimating the 
mass accretion rate assumes a negligible contribution from activity. This seems to be true 
for brown dwarfs on a statistical basis ('t he A8662 component occurs almost exclusively 
in accreting objects'; iMohanty et al.ll2005l ). and for 2M 1207, in particular, by the lack of 
measureable effects from chromospheric emission in other lines. 



According to lKoenigll (jl99ll ) the mass accretion rate of TTS is related to their magnetic 



M^/^jM:' ■ ■ rI' . We assume Rt ~ 2R^ for the disk truncation 
le et aPboOoh. The mass and radius of 2M 1207 are = 0.03 and 



field str ength, B 
radius (IMuzero 

i?* = O.27i?0 (jChabrier et al.ll2000l ). Taking account of the scaling factors, the average 
accretion rate of 2M1207 observed in May 2006 (M ~ lO~^°M0/yr) yields an approximate 
value for the surface field, B ^ 200 G. The same equation yields for a TTS with 0.8 Mq, 
1.5 R(7), and M ~ 10~^Mf7)/yr a s urface field of ~ 600 G. This validates the remark made 



by IScholz fc Jayawardhanal (120061 ) who, without citing numbers for the field strength, had 
argued that the magnetic field for a typical brown dwarf is expected to be roughly half 
that of a typical TTS. The magnetic field of TTS obtained this way is somewhat lower 



than measured values, which are typically in the kilogauss range (e.g. lYang et al.ll2005l . and 
references herein). There are, however, considerable uncertainties connected with the use of 
the Koenigl-relation, such as the simplified assumption of a dipolar field and the location 
of the inner disk radius with respect to the corotation radius that determines the scaling 
factor. Given these ambiguities, the order of magnitude agreement between the estimated 
and observed fields is plausible. 



4.3. Accretion versus activity in TWA brown dwarfs 



Given the presence of a magnetic field, inferred from the accretion rate, the absence 
of strong magnetic activity on 2M 1207 is remarkable. A very stringent upper limit to its 
X-ray emission was placed by a 50ksec Chandra observation (IGizis fc Bharatll2004l ). The 
radio emission was constraine d to the lowest value among the brown dwarfs in the TWA 
(jOsten &: Jayawardhanall2006l ) . In Table [3] we summarize the accretion and activity measures 
for all four TWA brown dwarfs. 



The data in Table [3] has been collected from the literature except for the X-ray data for 
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SSSPM1102. This object was identified as a probable su bstellar TWA mem ber in a recent 
proper motion study of tfie Supercosmos Sky Survey data (IScfiolz et al.ll2005bl ). SSSPM 1102 
is in tlie field-of-view of an arcfiived XMM-Newton observation (Obs-ID 0112880201). We 
ran the standard Science Analysis System source detection process, and find no X-ray source 
at the optical position of SSSPM 1102. To estimate an upper limit to its X-ray luminosity, 
we measured the counts at the expected position of SSSPM 1102 and in an annulus around 
this position (representing the background). After applying the appropriate area scaling 
factor and extracting the on-source time from the expo sure map, an uppe r limit to the net 
source count rate was estimate d using the algo rithm of iKraft et al.l (Il99ll ). The distance to 
SSSPM 1102 has been given bv lMamaiekl J2005h {d ~ 43 pc). This yields < 5.3 ■ 10^6 erg/s 
in the 0.5 — 10 keV band. 

It is constructive to compare 2M 1207 to TWA-5 B, the best-studied substellar object 
in the TWA after 2M 1207. Matching their properties shows the dichotomy of one brown 
dwarf that accretes but is not magnetically active, and another brown dwarf that is mag- 
netically active but does not accrete. Presumably, some mechanism suppresses actiyity in 



accreting substellar objects, as seems to be the case for classical TTS (jPreibisch et al.l 12005 



Telleschi et al.l 120061 ). Possible explanations that have been put forth in the context of clas- 
sical TTS are changes in the magnetic field structure and/or heating process of the corona 
as a result of accretion. The absence of magnetic activity on 2M 1207 is further surprising, 
as accreting substellar objects wit h similar magnetic fie l d strength, e.g. y 410Anonl3, have 
been shown to be X-ray emitters (IMuzeroUe et al.ll2000l : iGiidel et al.ll2006l ). For the substel- 
lar regime, the relation between accretion and X-ray luminosity has not yet been established 
on a statistical sample. If such a relation, as suggested by 2M 1207 and TWA-5 B, is con- 
firmed, this would provide one more piece of evidence for the TTS-like character of brown 
dwarfs. X-ray and radio studies of the remaining two TWA brown dwarfs, both of which are 
non-accreting, should be useful. 



5. Conclusions 

A substantial number of high-resolution spectroscopic observations has been collected 
for the brown dwarf 2M 1207. These data have allowed us to study its accretion variability 
on timescales from hours to ~ 3 yrs. While small variations are obviously present throughout 
the rotation (~ Id), major changes in the structure of accretion may take place over months 
or years. A summary of all Ha width measurements suggests a trend towards increasing 
accretion rate during the last three years. However, it is not clear if this trend reflects 
erratic changes or a systematic evolution. All in all, our dedicated monitoring has shown that 
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the accretion pattern on 2M 1207 is probably more complex than expected. It comprises, 
longterm changes of the accretion rate and (maybe independently) of the flow structure 
(months to years), variations on the timescale of the rotation rate (hours to days), and 
possibly on even shorter timescales as suggested by our high temporal resolution data. 

2M1107 is among the lowest mass dwarfs with measured mass accretion rate. The 
values derived from the various observations obtained in a time-span of ~ 3yrs, differ by 
nearly 2 dex, similar to (and possibly explaining) the scatter observed in the relation between 
accretion rate and stellar mass. Whether accretion variability is decisive for shaping the 
^acc ~ M relation will be established in future systematic spectroscopic studies across a 
wide mass range. 

This paper is based on data collected under the ESO program 077.C-0323. BS acknowl- 
edges financial support from ASl contract ASI-INAF 1/023/05/0. BS wishes to thank S. 
Randich and S. Mohanty for useful discussions on the UVES pipchne analysis and brown 
dwarf emission line spectra. RJ and AS acknowledge support from an NSERC grant to RJ. 

Facility: VLT:Kueyen (UVES) 
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Fig. 1. — Temporal evolution of parameters characterizing the Ha emission of 2M 1207 
during May 8-9, 2006. 
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Fig. 2. — top - Normalized average Ha line profiles of 2M 1207 obtained with UVES in May 
2006; bottom - Normalized variance line profiles and zero-variability level. 
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Fig. 3.— H« line profiles of 2M 1207 obtained with MIKE in Feb and Apr 2006. The data 
were normalized to the average continuum level near the line. 
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Table 1: Summary of Ha width measurements for 2M 1207. 



Date 




Reference 




n / 1 
[km/sj 




Apr 2002 


~ 300"f 


Gizis (^002) 


May 8, 2003 


170.. .200 


Mohantv et al. f2003) 


Jan 29 - Feb 1, 2005 


209. ..215 


Scholz et al. f2005a) 


Mar 17 - 19, 2005 


253.. .308 


Scholz et al. f2005a) 


Mar 27 - 30, 2005 


279.. .304 


Scholz et al. r2005a) 


Feb 21 - 22, 2006 


261. ..281 


present work 


Apr 11, 2006 


253 


present work 


May 8 - 9, 2006 


281. ..322 


present work 


May 16, 2006 


320 


Whelan et al. (2007) 



^ low-resolution data, no 10 % width measured. 
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Fig. 4. — Equivalent widths of Gaussian fits to the emission (WEQ;em) and absorption 
(M/EQ;abs) compoucut of the Ha hne, and normalized total observed flux (Fobs = — F^hs) 
versus Ha 10 % width. 
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Table 2: Equivalent widths and FWHM for detected emission lines. 



Identmcation 


^EQ 

r X 1 

[A] 


typ.error*- ' 
[A] 


[km/sj 


H/3 


4. 5. ..47.1 


~ 6 


180 


H7 


0.3. ..13. 5 


~ 2 


180 


m 






168 


He+CallH 






144 


CallK 






18 


HeIA5876 


1.8.. .9.1 


~ 2 


30 


Ha 


145.. .394 


~ 45 


210 


HeIA6678 


0.9... 2.4 


~ 2 


18 


CaIIA8498 


0...0.6 


~ 0.35 


12 


CaIIA8542 


0...0.3 


~ 0.25 


12 


CaIIA8662 


0...0.4 


~ 0.3 


12 



^"•^ No equivalent widths are given for the higher Balmer lines 
and for Ca II H+K because of the low continuum 
(b) Errors for the equivalent width are estimated 
by assuming a 15 % uncertainty in the 
continuum flux. 

^'^^ measured in an average spectrum that 
represents the typical line proflle for 2M 1207 
during the observation. 



Table 3: Activity and accretion measures of all known brown dwarfs in TWA. 



Designation 


IR excess 


W^Ha;10% 






FUV 


Refs 






[km/s] 


[erg/ s] 


[10^^ erg/s/Hz] 






TWA5B 




162 


41027 


< 2.0 




(-,1,2,3,-) 


2M 1207 




170.. .320 


< 1.210^6 


< 0.98 




(4,5,6,3,7) 


2M1139 




111 


? 


? 




(4,l,-,8,-) 


SSSPM1102 




194 


< 5.3 10^^ 


< 1.6 




(-,1,9,3,-) 



(1) - Mohantv et al. r2003) . (2) - Tsuboi et al. r2003) . (3) - Osten fc Javawardhana r2006) . 
(4) - Riaz et al. (2006), (5) - see references in this paper, (6) - Gizis fc Bharat (2004), 
(7) - Gizis et al. f2005l . (8) - Burgasser fc Putman f 2005^1 . (9) - this paper. 
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Fig. 5. — Mean of the standard deviation of fractional profile changes D{X) in the Ha hne 
as a function of time lag. Numbers on top of the panel indicate the number of data points 
in each of the time lag bins. 
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Fig. 6. — Correlation matrix of Ha for different subset of exposures (denoted i) obtained on 
May 8-9, 2006. left -t = l,..., 19; middle -i = 2Q, ...38; nght - i = 7, 10. 
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Fig. 7. — Contemporaneous Balmer emission line profiles (Ho; to He) for two epochs during 
the UVES run (8- frame averages). 
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Fig. 8. — Two spectra in the region of He I A5876 showing the change from the typical profile 
with extended red wing (solid line) to the weak, symmetric line (dotted line). 
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